Track reconstruction in the emulsion-lead target of the 
O ■ OPERA experiment using the ESS microscope 



o 

(N 
>^ 



(N 



>< 



February 1, 2008 



L. Arrabito^, C. Bozza^*^, S. Buontempo'', L. Consiglio'^, M. Cozzi^, N. D'Ambrosio^, 
^ ; G. De Lellis^ M. De Serio^l, F. Di Capua^ D. Di Ferdinando^, N. Di Marco^ 

^ '. A. Ereditato^, L. S. Esposito^, R. A. Fini\ G. Giacomelli^ , M. Giorgini^, G. Grella^°, 

; M. Ieva\ J. Janicsko Csathy*, F. Juget^, I. Kreslo^, I. Laktineh*^, K. Manai*^, 

G. Mandrioli^, A. Marotta^, P. Migliozzi^, P. Monacelli^, U. Moser^, M. T. Muciaccia\ 
Q ■ A. Pastore\ L. Patrizii^, Y. Petukhov"^, C. Pistillo^, M. Pozzato^, G. Romano^", 

, G. Rosa^, A. Russo^, N. Savvinov^, A. Schembri^, L. Scotto Lavina^, S. Simone^, 

>>■ M. Sioli^, C. Sirignano^o, G. Sirri^, P. Strolin^ V. Tioukov^ and T. Waelchli^. 

^ ! 

1. Dipartimento di Fisica deirUniversita di Bari and INFN, 70126 Bari, Italy 

2. University of Bern, CH-3012 Bern, Switzerland 

3. Dipartimento di Fisica dell'Universita di Bologna and INFN, 40127 Bologna, Italy 
^ ' 4. JINR - Joint Institute for Nuclear Research, 141980 Dubna, Russia 

04 ^ 5. Laboratori Nazionali del Gran Sasso, 67010 Assergi, L'Aquila, Italy 

O . 6. IPNL, IN2P3-CNRS and Universite Claude Bernard Lyon, 69622 Villeurbanne, France 

7. Dipartimento di Fisica deU'Universita "Federico 11" and INFN, 80125 Napoli, Italy 

8. University of Neuchatel, CH-2000 Neuchatel, Switzerland 

l/^ ■ 9. Dipartimento di Fisica deU'Universita "La Sapienza" and INFN, 00185 Roma, Italy 

' 10. Dipartimento di Fisica deU'Universita di Salerno and INFN, 84084 Fisciano, Salerno, Italy 

■ 

O ■ Abstract 



The OPERA experiment, designed to conclusively prove the existence of 
— » i^r oscillations in the atmospheric sector, makes use of a massive lead- 
5J] ' nuclear emulsion target to observe the appearance of Vr^s in the CNGS 

5^ I beam. The location and analysis of the neutrino interactions in quasi real- 

time required the development of fast computer-controlled microscopes able to 
reconstruct particle tracks with sub-micron precision and high efficiency at a 
speed of ~ 20cm^/h. This paper describes the performance in particle track 
reconstruction of the European Scanning System, a novel automatic microscope 
for the measurement of emulsion films developed for OPERA. 



1 Introduction 

In the last decades, several experiments using solar, atmospheric, reactor and accel- 
erator neutrinos have provided unambiguous evidence for neutrino mixing. Strong 

^Corresponding author. E-mail: Marilisa.Deserio@ba.infn.it 
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indications in favour of the hypothesis of Vr oscillations in the atmospheric 

sector have been obtained by [H [21 [3l 3] and have been more recently confirmed 
by E]. However, the direct observation of the appearance of oscillated z^^'s in 
the atmospheric signal region of the parameter space is still an open issue. The 
OPERA experiment [7j has been designed to achieve this goal using the CERN to 
Gran Sasso beam (CNGS [8]) over a baseline of 730 km. 

Based on the Emulsion Cloud Chamber technique [9] , nuclear emulsions [lOj are 
used as sub-micrometric space resolution trackers to detect the r leptons produced 
in Vt charged current interactions in a massive (0(kt)) lead target. 

The basic target unit (brick) with size 12.7 x 10.2 x 7.5 cm^ consists of 57 nuclear 
emulsion films, made of two 44 //m-thick layers on either side of a 205 fim plastic base, 
interleaved with 56 lead plates of 1 mm thickness. Bricks are arranged in vertical 
structures called walls. 

Electronic trackers, made of planes of plastic scintillator strips inserted in be- 
tween the target walls, and muon spectrometers accomplish the task of identifying 
in real-time the position of the neutrino interaction inside the target. As a result 
of the combined analysis of electronic data, a 3-D brick probability map is com- 
puted for each single event and the brick with the highest probability is promptly 
removed from the target during the run. The so-called Changeable Sheets (CS), 
pairs of emulsion films placed in front of each brick and acting as interfaces with 
the electronic trackers, allow the brick identification signal to be checked. Only in 
case of confirmation, the emulsion films of the brick are developed and sent to the 
scanning laboratories for interaction vertex location and further analysis. 

Assuming the nominal intensity of the CNGS beam (4.5 X 10^9 p.o.t. / year), 
about 30 bricks per day will be extracted from the target, corresponding to a total 
surface of nuclear emulsions to be scanned of a few thousands cm^ per daj0 In order 
to analyze the events in quasi real-time with a reasonable number of microscopes 
(~ 1 microscope / brick / day), the use of fast automatic systems with a scanning 
speed of about 20cm^/h is required. 

Moreover, the detection of the short-lived r leptons through the observation of 
their kink decay topology, namely an abrupt deviation of the r particle trajectory at 
the decay point in electron, muon and charged hadron(s), occurring over distances 
of ~ 1 mm at CNGS energy, demands adequate position (~ 1 /im) and angular 
(~ 1 mrad) resolutions. 

The European Scanning System (ESS) was designed to meet these experimental 
requirements The general structure of the system and the algorithms developed 
for image processing and particle tracking in single emulsion films, as well as the 
hardware and its performance, have been already described in detail [121 [T3| I14j . 
The performance in electron-pion separation has been recently published in [15\. 

In this paper, after a brief summary of the main features of the ESS, the film- 
by-film alignment algorithm, relevant for the analysis of OPERA bricks, and the 

^The position resolution of the electronic trackers 1 cm) will allow to reduce the area to be 
scanned to ~ 25 cm^ for CC interactions with penetrating muon tracks. For NC events, we foresee 
to scan the whole surface. 

different system for the scanning of OPERA films at high speed has been developed in Japan 
(SUTS [11]). 
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Figure 1: The European Scanning System microscope. 



tracking performance of the system are presented. The methods to be apphed for 
the location of neutrino events are then described and the results of a beam exposure 
to high-energy pions, designed to test the above procedures, are reported. 

2 The European Scanning System 

The ESS microscope, shown in Fig. [H consists of 

• computer driven horizontal and vertical stages equipped with high-speed pre- 
cision mechanics; 

• customized optics providing achromatic planar images of the emulsion; 

• high resolution and high frame-rate camera interfaced with a programmable 
frame-grabber and vision processor. 

By adjusting the focal plane of the objective lens through the whole emulsion 
thickness, a sequence of tomographic images of each field of view are taken at equally 
spaced depth levels, processed and analysed in order to recognise aligned clusters 
of dark pixels (grains) produced by charged particles along their trajectories. Each 
44 //m-thick OPERA emulsion layer is spanned by 15 tomographic images in steps 
of ~ 3 /im, accounting for the effective focal depth of the system. 

Micro-tracks are reconstructed as geometrical alignments of grains detected in 
different levels within the same layer. By connecting micro-tracks across the plastic 
support, the so-called base-tracks are formed (Fig. [2]). This strongly reduces the 
instrumental background due to fake combinatorial alignments, thus significantly 
improving the signal to noise ratio, and increases the precision of track angle recon- 
struction by minimising distortion effects. 

The reconstruction of particle tracks in OPERA bricks requires connecting base- 
tracks in several consecutive films. 
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Figure 2: Micro-track connection across the plastic base (base-track). 

3 Emulsion film alignment 

In order to define a global reference system (hereafter referred to as brick reference 
system), prior to track reconstruction a set affine transformations (shift, rotation 
and expansion) relating track coordinates in consecutive films have to be computed 
to account for relative misalignments and deformations. The mechanical accuracy 
of film piling in brick assembly is indeed of 50 -i- 100 //m, much worse than the 
achievable precision. Furthermore, emulsion films are affected by environmental 
conditions (temperature and humidity) altering their original geometry. The task 
of film alignment is accomplished by exposing each brick to a controlled flux of 
cosmic rays before disassembly, as described in [16j, and by applying the following 
procedure film by film. 

Starting from a set of independent measurements in single emulsion films, an 
aligned volume is created through an iterative pattern matching procedure comput- 
ing the parameters of the transformations 



where x^^^"^, y^*'™ are single film track coordinates and x^"'^^, y^"'^^ are the corre- 
sponding aligned ones. A least square fit is applied after maximising the number 
of matching pairs within predefined position and slope tolerances measured in three 
zones, typically chosen at the corners of the scanned area in order to maximise the 
lever arm and thus disentangle the contributions due to rotation and translation. 

Moreover, a wide angular spectrum of collected tracks allows the total thick- 
ness of the emulsions and of the interleaved lead plates to be precisely computed, 
accounting for possible a-planarities. The difference between the real position of 
a given track in the i-th film and the extrapolated position from the j-th film de- 
pends indeed on the uncertainty 5Z on the knowledge of the thickness and is a linear 
function of the track angle. By using a set of tracks crossing the brick at differ- 
ent inclination angles, a fit can be performed and the correction to be applied to 




the nominal value of the Z coordinate, perpendicular to the emulsion plane, can be 
estimated, thus improving the quality of the alignment. 

By applying the above procedure to each pair of consecutive films, relative dis- 
placements can be reduced to the level of a few fim, adequate for film-by-film track 
following (see Section 5), down to less than 1 ^m, as required by local-scale (~ 1 mm) 
analysis (this subject will be discussed in forthcoming papers). 

The achievable precision mainly depends on the number of penetrating tracks, 
proportional to the area of the measured zones. The density of passing-through 
tracks should be low enough in order not to spoil the topological and kinematical 
reconstruction of neutrino events; on the other hand, the scanning time is a critical 
issue and needs to be minimised. Typically, a density of the order of a few tracks 
/mm^ and scanning zones of 0(10 mm^) are a reasonable compromise between these 
two conflicting requirements. 

Once all films have been aligned, base-tracks are connected to form volume tracks. 
A two-step algorithm is applied [17] . The first step consists of identifying incremen- 
tal chains of consecutive base-tracks. These chains are then used as input for the 
Kalman filter algorithm [18], performing track fitting and propagation with a max- 
imum number of allowed consecutive holes (i.e. missing segments) accounting for 
tracking inefficiency. 

4 Study of the ESS basic performance 

In order to study the track reconstruction performance of the ESS in terms of res- 
olution and efficiency, a test exposure was performed at the CERN PS in the T7 
experimental area. One brick consisting of 64 emulsion films with no lead plates 
in between was exposed to lOGeV/c negatively-charged pions. The choice of high 
beam momentum and the use of emulsions only were motivated by the need to min- 
imise the effect of multiple Coulomb scattering that would spoil the measurement 
of the intrinsic resolution. 

The emulsion films, produced by the Fuji FilrrH company, were similar to those 
used for the experiment, although they were not part of the official OPERA produc- 
tion. They were transported by plane from Japan to CERN where the refreshing 
procedure was applied shortly before exposure. 

Nuclear emulsions are continuously sensitive to charged particles from production 
time until development. As a consequence, they normally accumulate a significant 
amount of background due to cosmic rays and ambient radioactivity that cannot 
be distinguished from particle tracks originating from the beam. For this reason, 
a process of controlled fading of the emulsions, called refreshing, consisting in the 
destruction of latent image centers in silver bromide crystals by keeping the films 
at high temperature and high relative humidity for a few days, was implemented 
for the OPERA experiment [19j. The working conditions can be tuned in order to 
reduce the average number of grains per base-track below the applied threshold, 
while preserving the emulsion sensitivity. This implies the cancellation of recorded 
tracks. The achieved track erasing efficiency can be higher than 98%. 

^Fuji Film, Minamiashigara, 250-0193, Japan. 
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angle (rad) transverse position resolution (micron) 

Figure 3: Longitudinal position residuals between base-tracks and corresponding 
volume tracks as a function of track angle (left); transverse position residuals (right). 



The angular dependence of the ESS tracking efficiency and precision were studied 
by rotating the brick with respect to the beam direction in 7 different positions from 
Omrad (tracks perpendicular to the films) to GOOmrad. The pion beam and the 
exposure time were adjusted so as to produce a density in emulsion of about 3 tracks 
/mm^ for each angle. Reference emulsions not exposed to the pion beam were also 
developed to provide an evaluation of the instrumental background, mainly due to 
uncorrelated grains randomly produced in the development process {fog) that tend 
to hide real track grains, hence resulting in fake coincidences. 

In order to collect a sample of several hundred tracks per angle exposure, an area 
of a few cm^ was scanned for each film and a pattern-matching procedure was applied 
to compute the parameters of film inter-calibration transformations, as described in 
Section 3. Measured base-tracks were then connected to form volume tracks. 

Fig. [3] shows the position residuals between each base-track and the correspond- 
ing fitted volume track as a function of the track angle. If S = (sx,Sy) is the slope 
vector of a straight line representing the projection of a particle track in the hor- 
izontal plane, the unit vectors ny = {ux^Uy) such that S = Sny and its normal 

= (n^, —rix) define two orthogonal directions, called longitudinal and transverse, 
that can be used to decouple the slope-dependent contribution from the intrinsic ac- 
curacy. The longitudinal component is indeed affected by a slope-dependent term 
proportional to the uncertainty on the Z positions of grains that can be assumed to 
be given by the sampling step of the tomographic sequence [12]. 

The left plot of Fig. [3]shows longitudinal residuals ranging from 0.3 /^m to 0.7 /^m. 
Transverse residuals (right plot) are independent of track angle and represent the 
ultimate precision (0.3 ^m) that can be achieved accounting for measurement errors. 

The angular residuals, computed as differences between base-track and volume- 
track angles, are shown in Fig. [H Perpendicular tracks (left) can be reconstructed 
with a resolution of about 1.7mrad; the resolution worsens with the angle up to 
6.3mrad for tracks crossing the brick at about 600mrad with respect to the perpen- 
dicular direction (right). 

In order to evaluate the tracking efficiency, minimizing the effect of base-tracks 
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longitudinal angular resolution (rad) longitudinal angular resolution (rad) 

Figure 4: Angular residuals between base-tracks and corresponding volume tracks 
for particles crossing the emulsions perpendicularly (left) and at 600 mrad inclination 
(right). 



not related to the beam, sets of six adjacent films were aligned and beam volume- 
track^ with at least five consecutive base-track segments were selected. Set by 
set, the efficiency was computed for the two external films as the ratio between the 
number of tracks with six measured segments and the number of tracks with five 
measured base-tracks only. By averaging over several films, the plot in Fig. [5] was 
obtained. The base-track finding efficiency is larger than 90% for perpendicular 
tracks; correspondingly, the micro-tracking efficiency, as given by its square roolU, 
is above 95%. 

The shape of the curve in Fig. [5] is related to the slope dependence of the average 
number of grains per base-track shown in Fig. [6j for small angles, the shadowing 
effect (see |12] ) enhances the signal, thus improving the efficiency. The effect falls 
off with increasing slope. However, for sufficiently large angles, the increase in the 
average number of grains due to longer path length in emulsion produces again an 
increase of the efficiency. 

The background measurement was performed using unexposed reference films 
refreshed and developed at the same time as the other films. A few hundred base- 
tracks selected applying standard cuts (see [12]) were visually inspected in order to 
separate the contribution due to residual cosmic rays from instrumental background, 
consisting of base-tracks with at least one micro-track generated by accidental com- 
binations of fog grains. About 1 fake base-track /cm^ was found in the angular 
range [0, 400] mrad. 

The measurements presented in this section were performed with the ESS work- 
ing at its maximum speed of ~ 20cm^/h. 

*Beam tracks were defined as those with a reconstructed angle within 3 a from the center of any 
of the 7 beam peaks. 

^Micro-tracks in the two emulsion layers of each film are independently recognised and the 
efficiency of the linking procedure can be assumed equal to 1 since it only depends on geometrical 
acceptances that can be properly tuned. 
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Figure 5: Base-track finding efficiency as a function of track angle. 
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Figure 6: Average number of grains per base-track as a function of track angle. 
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5 Particle track reconstruction in OPERA bricks 



Two different steps are foreseen in tlie event location strategy we plan to apply for 
the analysis of the OPERA bricks. The first procedure is the so-called scan-back: 
starting from a set of predictions provided by the electronic detectors, tracks of 
secondary particles coming from a neutrino interaction are searched for in the CS 
films; then, these tracks are followed back, film by film, from the most downstream 
emulsion of the brick to the interaction point where they originate. Whenever a 
track disappearance signal is detected (the track is not found in a certain number of 
consecutive films set according to the tracking efficiency), an area of about 25 mm^ 
is measured in several films around the candidate vertex point {volume scan) in 
order to fully reconstruct the event and apply topological and kinematical selection 
criteria. 

In order to test the system capability in reconstructing particle tracks in the 
OPERA emulsion- lead target, several bricks assembled with refreshed films were 
exposed to a low intensity 8GeV/c vr" beam at the CERN PS - T7 line. CS films, 
packed under vacuum in a separate envelope, were attached to the downstream face 
of each brick. 

About 1000 pions with inclination angles of about 50mrad with respect to the 
perpendicular direction were recorded in each target unit over an area of 9 x 9 cm^ . 
Before unpacking and development, CS films were removecj^l and each brick was 
exposed to cosmic rays for 6 hours in order to collect a set of passing-through tracks 
(density 1 track /mm^) to be used for film alignment. 

The OPERA event location procedure starting from CS films was mimicked us- 
ing the tracks found in both CS's and confirmed in the most downstream film of the 
brick as predictions. Film by film, each track was searched for in one single micro- 
scope view (390 x 310 /xm^) within a slope-dependent angular window, accounting 
for the decrease in accuracy with increasing angle. In case two or more base-tracks 
satisfied the angular cut, the closest-distance criterion was used for selectior0. Track 
propagation from film to film was performed by using the position and the slope of 
the most upstream track segment. Losses due to tracking inefficiency were recov- 
ered by requiring that the search for each track be iterated in 3 consecutive emulsion 
films if no candidate was found in one film. In case no base-track compatible with 
the prediction was found in 3 films, the position of the most upstream measured 
base-track was tagged as a candidate vertex point and the volume scan procedure 
was applied around it. 

The tracking efficiency, computed as the number of found base-tracks divided 
by the total number of films for the sample of passing-through particles (mainly 
primary pions and muons from beam contamination), is shown in Fig. [71 The 
measured value (~ 93%) is in agreement with the expected behaviour at small angle 
shown in Fig. O 

^Since the scanning of the CS films is a crucial step in the event location procedure (only tracks 
measured in the CS films are foUowed-back across the brick), the background track density should 
be kept as low as possible. 

^Different best-candidate selection criteria (for example, criteria based on the angular agreement) 
are currently being tested. 
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Figure 7: Tracking efficiency for non-interacting particles. 
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Figure 8: Left: Number of base-tracks per view. Right: Number of base-tracks per 
prediction satisfying the selection criteria. 



On average, 2 base-tracks with space angles smaller than 400mrad were recon- 
structed in each microscope view around the position of a given predicted track (Fig. 
[8l left). By requiring that the slope difference between found and predicted base- 
tracks be not larger than (0.03 -|- 0.05 x slope), hence allowing for the non-negligible 
contribution due to the scattering of low energy particles in lead, the average number 
of candidates per prediction reduces to 1, the fraction of predictions with multiple 
candidates being of the order of 10~^ (Fig. El right). 

We foresee to complete the scanning of each brick in a few hours. In order 
to minimise the scanning time required for event location, fine corrections are not 
applied in scan-back mode. First of all, only a global film inter-calibration procedure 
is performed. Fig. [9] shows the distribution of the position differences between 
predicted and found base-track segments. The plot refers to tracks scattered over 
an area of about 40 cm^. The sigma of the gaussian fit is of about 7/im. Moreover, 
only one microscope view is measured for each predicted track, thus effects such 
as local planarity cannot be corrected. The resulting angular precision is of about 
4-^5 mrad. 
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Figure 9: Position differences between predicted and found base-track segments. 
Tracks are distributed over an area of about 40 cm'^. 
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Figure 10: Position (left) and angle (right) differences between base-tracks selected 
in two scan-back procedures corresponding to the same set of predictions. 



The scan-back reproducibility was studied by performing the procedure twice for 
the same set of predictions and by comparing film by film the selected base-tracks 
per each predicted track. As shown in Fig. \T0\ the reproducibility is of about 2.6 fiui 
in position and 3mrad in angle. 

As a cross-check for the validation of scan-back results, the volume scan pro- 
cedure was applied around track disappearance points and, by a combined analy- 
sis of scan-back and volume scan data, events were classified as single and multi- 
prong interactions, low energy (scattering) tracks and passing-through tracks. Film- 
dependent effects, such as local distortion and misalignments resulting from large 
rotations, not completely corrected in the inter-calibration procedure, were found 
to be the main sources of inefficiency (< 5%) in the track following. Such effects 
can be accounted for by enlarging the search window in case no base-tracks com- 
patible with a given prediction are found and by implementing a prediction-forking 
method, allowing for multiple candidate following, if the slope/position acceptances 
are increased. Failures of the system in passing-through track following can be thus 
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reduced to ~ 1%. 



6 Conclusions 

The particle tracking performance of the European Scanning System (ESS), a novel 
high-speed automatic microscope for the measurement of nuclear emulsion films, has 
been presented. 

The system, developed for the location and analysis of neutrino interactions in 
the emulsion-lead target of the OPERA experiment, fully satisfies the requirements 
of high tracking efficiency and accuracy, working at a speed of ~ 20cm^/h. 

The results of a test exposure to high-energy pions, designed to study the meth- 
ods to be applied for film-by-film alignment, track following and event location and 
confirmation, were discussed in details, showing that the system is able to track par- 
ticles produced in beam interactions with lead through the whole brick and identify 
the vertex points where they originate. The achieved results make us confident in 
view of forthcoming CNGS neutrino runs. 
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